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Abstract 



We explore the duality between supersymmetric Wilson loop on null polygonal contours in 



maximally supersymmetric Yang-Mills theory and next-to-maximal helicity violating (NMHV) 
scattering amplitudes. Earlier analyses demonstrated that the use of a dimensional regulator 
for ultraviolet divergences, induced due to presence of the cusps on the loop, yields anomalies 
that break both conformal symmetry and supersymmetry. At one-loop order, these are present 
only in Grassmann components localized in the vicinity of a single cusp and result in a universal 
function for any number of sites of the polygon that can be subtracted away in a systematic 
manner leaving a well-defined supersymmetric remainder dual to corresponding components of 
the superamplitude. The question remains though whether components which were free from the 
aforementioned supersymmetric anomaly at leading order of perturbation theory remain so once 
computed at higher orders. Presently we verify this fact by calculating a particular component 
of the null polygonal super Wilson loop at two loops restricting the contour kinematics to a 
two-dimensional subspace. This allows one to perform all computations in a concise analytical 
form and trace the pattern of cancellations between individual Feynman graphs in a transparent 
fashion. As a consequence of our consideration we obtain a dual conformally invariant result for 
the remainder function in agreement with one-loop NMHV amplitudes. 



The duality between scattering amplitudes in M = 4 super- Yang-Mills theory and a super- 
symmetric extension of the Wilson loop spanned on a polygonal closed contour with its sites 
defined by particles' momenta involved in scattering occupied an important niche in devising 
new techniques for analysis of dynamics of gauge theories at weak and strong coupling regimes 
and interpolation between the two. A distinguished role played by the maximal supersymmetry 
in four dimensions is that all particles of the theory can be combined into a single CPT self- 
conjugated light-cone superfield $ defined by a (finite) series in the Grassmann variable r\ A with 
coefficients determined by the fields of appropriate helicity to compensate for the deficit intro- 
duced by the r\ itself and matching SU(4) tensor structure [U12]- Thus, the n— particle S-matrix 
of the theory is concisely represented by the amputated Green A n functions of n superfields $. 
Extracting the (super) momentum conservation laws allows one to cast the superamplitude A n 
into the following form [2] 

An = * (2?r) (12) (23) (34)... (n- In) MK K Vi) " (1) 

The use of the spinor helicity formalism, adopted here and below, simplifies the representation 
of the amplitude. Namely, the massless particles' momenta pf a = X a X a and their chiral charges 
gA,a _ yjAya are wr j^ en by means of commuting Weyl spinors X a and A", with their inner 
product defined as (ij) = \f\j a and [ij] = A i(i A", as well as anticommuting Grassmann variables 

r/ A transforming in the fundamental of SU(A). The reduced amplitude A n admits an expansion 
in terms of r/'s 

A n = A n ,0 + A n ,l + ■■■ , (2) 

that terminates at ordei0 k = n — 4, with each term being a homogeneous polynomial of degree 
f] 4k . Each term in this expansion describes scattering of particle with total helicity — n + 4 + 2k, 
with the leading term being the maximal helicity- violating amplitudes (MHV), then the next-to- 
maximal helicity- violating (NMHV) amplitude etc. At tree level, A n °\ = 1, while the latter can 
be written as a sum [3] 

A n \ — Rn;qr > (3) 

Kq<r<n 

of superconformal invariants 

R = £ 4 ((n, q-l,q,r- l)xr + cyclic) 

(q — 1, q, r — 1, r){q, r — 1, r, n){r — 1, r, n, q — n, q — 1, q)(n, q — 1, g, r — 1) ' 

written in terms of momentum twistors Z'- = (A", a;" Q Aj Q ), with angle-brackets being (ijkl) = 

e a bcdZ?ZjZ%.Zf. In what follows, the focus of our analysis will be A n ,i at higher orders of per- 
turbation theory. 

A profound realization of the past few years was that the superamplitude ([T]) is expected to 
admit a dual representation in terms of a Wilson superloop [H E] spanned on a closed polygonal 
supercontour with its vertices localized at (xj, 9f) such that its path segments are proportional 
to the particles' (super)momenta, pf a = (xi — Xi + \) aa = xff +1 and q Aa = (9i — 9i + i) A,a = ii !^ 1 , 

(W^A)) = ^(tr(W [ln] ...W [32 ]W [21] )) , (5) 



^^Nilpotcncc of the Grassmann variables alone is not sufficient to produce this constraint, the reduction by 
maximal degree by four is a consequence of superconformal symmetry. 
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where 




(6) 



V Jo / 

where the path is parametrized by X[u + i](t) = Xi — txu + i, 9[u + i](t) = 9i — t9a + \ and the first few 
terms in the superconnection read5 



It is a generalization of the duality for the lowest k = MHV component 0, El El EE] elucidated 
by now through multi-loop calculations [HI [121 1131 EH H3 EE E] • The above equation has the 
unusual property of mixing orders of perturbation theory on both sides of the equation, for 
instance, the £-th order k = 1 NMHV amplitude emerges from a (£ + l)-loop computation of the 
superloop, however, from terms quartic in Grassmann variables, etc. 

To lowest order in coupling the amplitudes and their duals on the super Wilson loop side are 
expected to be invariant under the so-called dual superconformal symmetry which acts on the 
dual coordinates (xi,9i) [3]. At sub leading orders in coupling, some of the symmetry generators 
are broken by the ultraviolet regulator in a predictable fashion. This was clearly demonstrated 
in great details for MHV amplitudes and its dual bosonic Wilson loop in Ref. [9]. However, 
the first encounter with the superloop's rj 4 — component, dual to the tree NMHV amplitudes, 
demonstrated that the former suffer from another anomalous effect [18]. Namely, the use of the 
Four-Dimensional Helicity scheme [19], adopted for the bulk of higher loop calculations on the 
amplitudes side as it preserves the spinor-helicity formalism, induces a conformal and super- 
symmetric anomaly which breaks the above correspondence [TBI [20] • However, this anomalous 
contribution has a universal form and can be subtracted away in a consistent manner, restoring 
the supersymmetry and conformal symmetry and thus resuscitating the conjectured duality. It 
is important to realize at the NMHV level, the degree four Grassmannian structure becomes 
anomalous provided it contains at most three adjacent indices [20], e.g., Xi-iXiXi+ij Xi-iXi e ^ c - 
Therefore, any structure where at least one of the indices is not adjacent to the rest will be 
conformal and given by the corresponding component of the /^-invariants. The question still re- 
mains whether those components that were not anomalous at leading order develop unexpected 
anomalies once computed at subleading orders. This is the issue that we will address in the 
present study. 

We will perform a two-loop computation of a non-anomalous component at leading order, 
picking X2X3X6X7 as the object of analysis. In order to be able to track explicitly all intricacies 
of cancellations between Feynman diagrams without the complications of dealing with higher- 
degree transcendental functions intrinsic to computations in the full four-dimensional kinematics, 
we will restrict the contour of the superloop to a two-dimensional subspace [61 [211 122] . In this 
situation, the highest transcendentality that one can expect in the result is degree two, which 
encompasses dilogarithms and squares of logs (as well as their lower powers). Moreover, the 

2 Thcsc arc the only components that we will need for the main calculation performed in the paper. 




The duality relation between the two objects is of the following form 
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Figure 1: Null octagonal Wilson loop contour in the two-dimensional kinematics. The exchanged 
scalar between the cusps x 3 and £7 selects the one-loop X2X3X6X7 component of the superloop 
and expressible in terms of the corresponding Grassmann projection of the -Rs ; 37 superinvariant. 



first nontrivial loop has the octagonal shape as shown in Fig. [TJ along with the definition of the 
light-like directions of the segments. The component in question of the tree NMHV amplitude 
is expressed in terms of i?8;37, 

R X2X3X6X7 /qx 

^ 37 " 2(23) (67)^3 [ ) 

This results can be easily reproduced by evaluating the corresponding component of the one-loop 
super Wilson loop, which is determined by the correlation function of two superconnections B 
each of which gets reduced to boundary terma^l with the scalar field localized at the vertices x 3 
and x 7 , 

/VA,(1)\ 9 2 C F X2X3X6X7 , 2 + - \£ / 1n \ 

( = "4^ (23X67)^3 ( ^ ^ ' (10) 
Here we kept the regularized form of the one-loop result since it will be essential for the definition 
of the remainder function in the discussion that follows. Notice that we absorbed transcendental 
constants into the rescaled mass parameter 2ne yE fi 2 — > fi 2 . Removing the regulator, e — > 0, we 
immediately see that the one-loop superloop is expressible in terms of the i?8;37 component of the 
superconformal invariant. This is the expected result since the anomaly emerges only in adjacent 
Grassmann components as explained above. 

The complexity level of the computation that follows is comparable to the two-loop calculation 
of the bosonic Wilson loop which is dual to MHV amplitudes. Presently, the two-loop analysis 
yields the dual to the one-loop NMHV amplitude since we are extracting degree-four Grassmann 
component. As in our previous studies [T5j |2"U] we will adopt the Four-Dimensional Helicity 
scheme [19] to regularize divergences in Feynman graphs. This regularization is the closest one 
to the way one tackles infrared divergent scattering amplitudes. The details of the analysis are 
deferred to the Appendix. 

Due to the choice of the particular Grassmann component, a number of Feynman graphs 
should not taken into account. Namely, at second order in coupling, one has to include the 



3 For instance, W[32] = W[4 3 ] = 1 + \g 2 X2Xz 0ab(£3)/(23) and analogously for other two segments adjacent 
to the vertex x-j. 



3 




Figure 2: Cancellation mechanism between the different ordering of gluon and scalar emission 
with the seagull terms stemming from the field component of the covariant derivative D4> A b- 

effects from the covariant derivative (see the last term in Eq. (J7J)) along with emission of scalar 
and gluon fields off the super- Wilson links. However, a quick inspection demonstrates that the 
sum of two orderings of emission along with seagull terms vanish as shown in Fig. [2j The 
cancellation works as follows. Consider the [32]-superlink as an example. Expanding it to second 
order in g, and keeping track of X2X3 component only (and ignoring fermions for a moment), we 
find 

Vita X = ^^M-l f 1 dtt(2\[A(t),4> AB (tm\ (11) 



4 (23) 



'0 



1 nt 

dt / dt' 



Jo 



(2\A(t)\2\ [t'± + lj 4>ab (0 + (tj t + lj 4> AB (t)(2\A(t')\2] 



Here the argument of all functions involved stands for f(t) = f(x[2s\(t)). The first line above 
displays the gauge field part of the covariant derivative, while the terms involving derivatives in 
the second line emerge from its flat part. Finally, derivative- free contributions in the integrand 
of the two-fold integrals come from two ordering of inserting the gluon and the scalar field into 
the [32]-link. The follow-up simplification of this expression is straightforward and one finds that 
the scalar fields is nailed down to the vertex at x = 23 while the gluon is emitted from any point 
on the link 

W [32] X ^ -^^4> AB {x,) j\t{2\A{t)\2\. (12) 

It takes the form of the leading order scalar emission vertex and a bosonic Wilson segment 
attached to it. Analogous arguments apply with minor modifications to other superlinks adjacent 
to the cusps at x 3 and X7 yielding contributions with the scalar localized at the cusps and gluon 
strings attached to it. 

As a consequence of this consideration, we are left with graphs where the scalar can only spill 
off the cusps at 23 and X?, thus generating two- loop diagrams shown in Fig. [3] (a-g) along with 
other attachment of gluons to other segments of the contour to form a gauge- invariant set. Since 
the same Grassmann structures, either X2X3 or X6X7> can t» e induced by fermions emitted off the 
two adjacent links, to the the given order in coupling, there is an extra graph of the type (h). 
Notice that it is required by supersymmetry of the superloop and it will be instrumental for the 
cancellation of the double-pole divergences in the non-abelian color structure CfC a . 

To present the result of our analysis, we will strip the dependence on the Grassmann variables 
and powers of the gauge coupling constant from the component of the super Wilson loop that 
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we are interested in 



4tt 2 y (23) (67)x+ 7 x 



37 „ V 



"IK 



-C F C A w? 



(13) 



Here the sum runs over the diagrams displayed in Fig. [3] and split it into abelian and maximally 
non-abelian color Casimirs. The contribution to the abelian part of the Wilson loop stems from 
the diagrams in (a), (b), (c), (d) and (e) and the result of rather elementary computations gives 



w 



(a) 



2 ( X 73' E 73/ i ) [( x n X &Z^ ) + ( a; 72 :r 87/ i ) 



,2\£ 



xi 6 x u ^r (14) 

X 36 ;r 34/ i2 ) £ + ( — x 23 X 83f J ' 2 y] { £ ^ + C2) 7 



W 



(c) 



W 



(b) 



W 



I In 


^83 


ln^f 


-Ii„ 


2-74 


In ^Z3 
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x 73 


x 72 
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^73 


1J.1 _l_ , 

x 63 


NA 






i/i NA 
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(b) > 




m k = 


«>(d) , 




w (e) = W 



NA 

(e) 



(15) 
(16) 



with w^ d e -) displayed below. The abelian part of the expression has the following multiplicative 
structure 2 



X2X3X6X7 

(23)(67)x+x 



37 



where first factor is the one-loop correction to the superloop (Wg 1 ^) defining the tree NMHV am- 
plitude and the second is the one-loop correction to the bosonic loop which is equal to Cp J2 a w a 
up to factors of the coupling constant. Thus the remainder function, defined by subtracting the 
ultraviolet divergent contributions conventionally by 



^l = (wg>-<wg><wg> 



;(2) 



(1) 



is solely defined by the maximally nonabelian color. Therefore, the sum of all corresponding 
contributions has to be dual conformally invariant. 

The analysis of the maximally non-abelian contributions is more involved. The result of 
rather lengthy calculations can be cast to the following form 



w 



(b) 



1 

2< 



-x+x 73 /i 2 ) £ [(-x+x^Y + (-x+a; 87 /i 2 ) £ ] (e 2 + C2) 



w 



(d) 



1 

2 



x 



62 



-78 



73 



X 



71 



^ In ^ In ^ + i In ^6 in ^1 + i In ^63 in % + i In ^23 in 



x 



62 



X 



31 



73, 



;r 



(19) 
(20) 



83 



W 



(e) 



■ft' 



(f) 



iln^ln^ + -ln X 

•^67 "^48 



i>2 j ^34 



23 



X 



8-1 



£ ( Xj 3 X 73 /l ) 



+ 



"1 


j x 73 


ln^ 2 +1 

x 73 2 


X 73 


2 


2x 23 _ 


zx 76. 



ln^f 



1 + 



73 



2x 



87 



In ^3 



+ 



x 



73 



1 + 



73 



X 



34 



In^Zi 



73 

+ 2 



73 



(21) 
(22) 
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+ 



+ 



X 



73 



2^23 



x 



"2.3 

+ 
73 



+ 



X 



73 



2^76 



X 



X 



X 



73 



2x 



Li 2 



x 
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2x 



23 



ln 2^72 



X+ + 4 



X 



73 



2x 



76 



34 

In 2 ^ 



x 



x 



X 



73 



2x 



Li 9 



87 



x 



- 2 



x 



73 



1 

+ 4 



1 + 



x 



73 



2x 



87 



X 



x 



73 



1 + 



X 



73 



2x 



34 



ln 2 ^i, 



x 



w 



(g) 



-x 72 x 73 /i 



2^ x _73 j n x 72 



X 



23 



X 



+ + a \ X &3 X 73^ 
73 ^ 



73 



2^£ X 73 i ^63 
+ + 



(23) 
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l(-xf 3 x 8 V)^ln^ - i(-x+x^ 2 )^ln^ 



x 



87 



"73 



X 



31 



X 



73 
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31 
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+ 



Sn^34ln^7+ln^8ln^ + ln^ln X 



67 
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+ 
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X 



73 



2x 



23 



Li 9 
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2x 
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x 
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X 



73 



X 
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73 



+ 



1 + 
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73 



2x 



34 



x 
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(h) 



2\e 



+ 2~( a '73 ;r 73/ i2 ) e ( x 76 X 87f j2 ) 6 



e- 2 + C 2 



<T 2 + C2 



x 72 x 74 

X 73 X 7 3 



(24) 



r -l ln W83 
X 73 a; 73 J 
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4 



- j n 2 X 72 X 7Z 



1 

4 



2 J R 2 • x '63- x '73 



< = -e" 1 ( 



73- 6 74 x 73 x &3 
x '73 X 73f J ' ) — ^ . 
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Li 2 l^ 



x 



73 



X 
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x 



73 
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Li 2 (^ 



(25) 



Summing up the diagrams, we find that, as anticipated, that all ultraviolet poles chancel 
among different graphs in the remainder function. In particular, the fermionic graph (h) exactly 
cancels the maximally non-abelian color structure in diagram (b). The sum of single poles vanish 
as well provided we set (—xfjX^fi 2 ) 6 — > 1, but otherwise they induce finite contributions. The 
latter are of paramount importance in removing all terms proportional to squares of the logs 
that depend solely either on plus or minus components of the dual coordinates. Finally, all 
dilogarithms present in individual graphs sum-up to zero in the remainder function as well. As 
a result the latter takes a very simple form 

1 ( ( ' ■ X2X3X6X7 

(26) 



*s = ( 



9_ 

4rr 2 



(23)(67)x+x 73 

x [\nu + \nu~ + \nu + \n(l + u~) + \n(l + u + )\nu~ -\n(l + u + )\n(l + u~) + 2( 2 ] , 
upon the introduction of the conformal cross-ratios 

• L 32 X Q7 - x 87 x 34 



U = 



U 



x m x 73 



X 84 X 73 



(27) 
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(a) (b) (c) 




(g) (h) (i) 



Figure 3: All topologies of Feynman diagrams contributing to the X2X3X6X7 component of the 
supersymmertic Wilson loop at two loop order. The blob on the scalar line in (i) stands for the 
sum of vacuum polarization bubbles due to gauge fields and gauginos. 

This expression is an agreement with the result of a recent analysis that bypasses the calculation of 
the Feynman graphs and finds the result in question by integrating the Ward identities associated 
with Q supersymmetry from the the tree- level NNMHV amplitude [23J. 

Presently, we verified by a brute-force Feynman graph calculation the duality between the 
supersymmetric extension of the null polygonal Wilson loop and the super amplitude in maximally 
supersymmetric gauge theory. Our analysis elucidates the validity of the correspondence for 
Grassmann components which do not involve at least three adjacent particle indices. The latter 
were shown to be anomalous already at one-loop order. However, once the universal conformal 
anomaly is subtracted out, the duality gets restored. In the forthcoming work [24j, we will 
demonstrate how one can perform a super-gauge transformation on the super Wilson loop in 
order to gauge away in a systematic manner the notorious anomalous contributions. 

We would like to thank Gregory Korchemsky for usefull correspondence and instructive com- 
ments and Song He for making the component form of NMHV superamplitudes in restricted 
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kinematics available to us. This work was supported by the U.S. National Science Foundation 
under grants PHY-0757394 and PHY-1068286. 



A Details of the calculation 

The bulk of graphs is easy to compute. Let us pay special attention to a couple of them that are 
not as straightforward, namely diagrams (g) and (h) in Fig. [3] 

A.l Diagram (g) 

Using the usual Feynman rules, we find the following integral representation for the graph (g), 

W 8;l W = » (4^)3 *12 (23)(67) r (1-^(2-0 J o dtJ g (t), 

where 

[z-x-jY (z-x 3 y 



Jg(t) = / d^ 2e Z 



[z — X7) 2 (z — x 3 ) 2 



\ z - X[12] ( t ))Y- £ [-(z - x 3 n-c[-( Z - X7 n- 

(28) 

As a first step, we use the Feynman parametrization to put J g into the form 
Jit) = - ^(2-2.) (29) 



X /i 



,.. Z" 1 dsi ds 2 ds 3 5 (YS=i s i ~ !) (sis 2 s 3 ) £ r + + + 

/ 7 2 2 2 i2-2e S l S 3^[i2] 3 _ S l S 2^[i 2 ] 7 + 2S 2 S 3 X 73 

JO [ — Si^^p^n — S l' S 3 a; [i2]3 _ s 2'S3X 7 3j 



Now, expressing S3 in terms of the other two variables via the S— function constraint and changing 
the integration variables aa^ s 2 — > Sis 2 and as a consequence s 3 — > sis 2 , we can integrate with 
respect to t to get 

1 iir 2 ~ 6 T(l - 2e) (-2fi 2 ) 2e f 1 ds\ f 1 ds 2 sis 2 x% 3 - s x s 2 x~l 7 + 2s x s 2 s 2 x^ 3 

a{) ~ Y 2 {l-e)Y{2-e) Ax~ 12 J, ^ J (s^F s 2 xf 7 + * 2 s+ 

1 1 1 



[sis 2 x 71 x 71 + sis 2 x 12 x 23 + sis 2 s 2 x 73 x 13 } 1 2e s\ 2£ {x 73 Y~ 2£ [s\x 71 + sis 2 x 73 Y 2e 

(30) 

It might appear that there is a double pole emerging from this integrals, however, the only pole 
that realizes is from the s 2 integration in the vicinity of s 2 = 0, while the expression in curly 
brackets scales as s x and it tends to zero canceling the potential singular behavior. We extract 
the pole in the s 2 integral via the following formula 

1 tl s _i ,/ n s f 1 ds 2 



^- e f{s 2 ) = s-V(O) + / T^LfM - /(0)] . (31) 

s 2 JO s 2 



^ere and below s = 1 — s for any variables. 
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This allows us to cast the result after some manipulations into the form 



dtJ g (t) 



i7r 2-eT(l-2e) 



1 



with the set of ij integrals that can be evaluated with the result 
^ dS! 1 



[e-\-2^x 73 ) 2e h + ~ 2sn) J 2 + x+ 3 h] , (32) 
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73 
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73/ J 



(34) 
(35) 



Notice that the second integral involves a denominator that mixes both plus and minus com- 
ponents. This effects disappears one we add a mirror symmetric diagram yielding a factorized 
product of function of plus and minus variables. Summing all diagrams of this topology we get 
the expression in Eq. (123]) . 

A. 2 Diagram (h) 

Now, we turn to the second graph. 



W (2) 

yV 8;l(h) 



- ^ 4 ^_^ X2X3X6X7-|p|-r(l - g)r 2 (2 - e) J ds dtJ h (t,s) 



where we have used the identity [2\(z — x 2 z(t))(z — X[ 34 ](s))|3] = [23] (z — X3) 2 in order to define 
the coordinate integral 



Jh(t,s) = fl 



it- 



d A - 2e z-, 



(z - x 3 f 



(z - x 7 fY-*[-{z - x m (t))^[-(z - x m {s)) 



212-E ' 



(36) 



By means of the standard Feynman parametrization, one can cast it in the form after integration 
over z 

m 2 - £ Y{2 - 2e) 



where 
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h{t, s) 
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(37) 

(38) 
(39) 



To perform the integrations efficiently, we remove the S3 variable with the 5-function and then 
rescale S2 — > sis 2 which implies S3 — > sis 2 . The denominator admits a factorized form with two 
factors both linear in s%, Sis 2 x 2 ^ + Sis 3 x 2 ^ + S2S 3 xp 23 j j 34 j = 2s\ [Asi + B] with 

^ = (^73 - ^2^23) ( X 73 + s ^2 x m) , B = S 2 S 2 tS X 23 X 34 ■ ( 4 0) 

The next integration to be performed is with respect to s±, which produces a hypergeometric 
functions 2-^1- However, we notice that the small-Si region yields a contribution inverse in B 
such that the two subsequent s- and t-integration induce divergencies. This allows us to extract 
the leading inverse-power behavior of the regularized integral at small B and then resum the rest 
for e = 0. This gives 

( l te lS \-° _ r(i- £ )r(2- e ) 1 1 1 

J [A Sl + B}^ r(3-2e) B l ~ e A 2 ~ e 2A{A + B) 2 2A 2 {A + B) y) ' { } 

Then integrating over the s and t variables, we get the following representation for the integral 

■1 rl 1 rl rl ,1 f 1 1 1 

(42) 



f dt f dsli(t,s) = — I dt j ds f ds 2 s 2 s 2 
Jo Jo 16 J J J q 



_A(A + B) 2 A 2 (A + B) 
2fi 2 xt 3 x^r(-2fi 2 x+ 3 x^r T(l - e)F(2 - e) f 1 



8 x 23 a '34( a; 73 ;E 73) 2 ^(3 2e) 

where the divergent one- dimensional contributions are given by 



ds 2 S £ (s 2 )T £ (s 2 ) , 



*:^)= I 4 M . -= ~- r^r , Us 2 ) = f *__L__. (43) 



s 



' 6 [1 + SS 2 X 34 /x 73 } 2 6 J t 1 6 [1 — ts 2 X 23 /Xj 3 } 2 



Their e— expansion is easy to construct and reads to order 0(e), e.g., for T e (s 2 ) 
T e {s 2 ) = \ + - In (l - *S) (44) 



£ X 73 X 23 S 2 \ X 73 



+ e 



X 23 b 2 Z,X J3 X 23 b 2 I X 23 \ 1 2 / X 23 \ / X 23 



In l-s 2 ^ --lnMl-s 2 ^ -2Li 2 s 2 

. x 73 X 23 S 2 ^73 ^S 52 V x 73 / ^ V x 73 / V x 73 J . 



and the one for S £ (s 2 ) being analogous with the obvious substitutions of the defining variables. 
With poles being extracted explicitly, the remaining integrations can be performed with Math- 
ematica. The output is given however, in a form that involves dilogarithms with arguments 
depending of products of plus and minus variables. Instead of relying on known identities be- 
tween the dilogarithms to simplify the result and cast it as sum of functions depending either on 
plus or minus variables, the use of the formalism of symbols fZE\ ITS] becomes very instrumental 
for fast and efficient derivations of the sought identities. Just to give an example, we encounter 
the following combination of dilogarithms in the output, 

L(u,v) = Li 2 (l + — ) -Li 2 1 + — -U 2 (uv) , (45) 
V uv/ V u J 

where u = x 23 /x 73 and v = x± 3 /x-~ 3 . In order to disentangle the u and v dependence, we calculate 
the symbol of the right-hand side of this identity and find after simple manupulations 

V f V \ VU f vu \ 

S \L(u, v)} = = ® 1 1 + — J H ® H ) + (u + vu)® (uv) (46) 

uv \ uv/ u V u / 
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= (u + vu) £g> (uv) + (uv) (g) (m + f w) — t> eg) u — u®v — V ®v — U ® U + V ® V . 

From here, we can immediately read off the expression for the function itself (with a poten- 
tially present additive transcendental constant fixed by comparing both sides of the equation 
numerically) , 



L(u, v) = \n(uv) ln(w + vu) — In v In u 



In 2 v + Li 2 (w) 



Li 2 (v) . 



(47) 



The same technique is applicable to all other terms. The sum of all terms yields expressions with 
arguments being functions of either u or v variables separately. 
Finally, we get for the integral ii, 



dt / ds I\ (t, s) 



-2/x 2 x+3X3 4 ) £ (-2/i 2 x+x7 3 r r(l - e)T(2 - e) 



' x 23 x 34l x 73 x 73J 



T(3 - 2s) 



J__Il n ^2 _I ln ^Z4 



e 



x 



73 



X 



73 



+ 
+ 



ln{xj 2 /xj 3 ) 



ln(x 74 /x 73 ) 



' x 23 x 34l x 73 x 73J 



16 (^23) 2a '34 a; 73( a; 73) 2 ^ x 23( x 34) 2 ( x 73) 2x 



73 



In (^72^73/^73^34) ^(a^A-re) 



JZ x 23 x 34l x 73 x 73J 



' x 23 x 34l x 73 x 73/ 1 



U 2 (x i3 /x 73 ) 

' np i~y> ( np np ~\ 2 

' x 23 x 34l x 73 x 73J 



(48) 



The calculation of the second contribution I 2 is much simpler since all one is after is the 
double and single pole part of the integral since they get compensated by the overall factor of 
e. Then in the analogous to Eq. (|4ip integral with respect to Si, one has to keep only the first 
term. The subsequent integrations over s and t like done above in Eq. (144]) . immediately gives 
the final answer 



dt 



'dsu^s) = t3^t3x u n-2f, 2 xj 3 xj 3 r r 2 (i- £ ) 



A np'' np np~^~ np 

% x 23 x 34 x 73 x 73 



T(2 - 2e) 



x 72 j n ^72 



+ 



X 



74 j ^74 



X 



23 



x 



73 



x 



13 



^73 



(49) 

Summing both contributions together, we find half of the result displayed in Eq. (I24p . The other 
half is given by the mirror symmetric diagram, computed via the formalism outlined above. 
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